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ABSTRACT 

Recent observations of z ~ 6 quasars and the cosmic microwave background imply a complex history to cosmic 
reionization. Such a history requires some form of feedback to extend reionization over a long time interval, but 
the nature of the feedback and how rapidly it operates remain highly uncertain. Here we focus on one aspect of 
this complexity: which physical processes can cause the global ionized fraction to evolve non-monotonically 
with cosmic time? We consider a range of mechanisms and conclude that double reionization is much less likely 
than a long, but still monotonic, ionization history. We first examine how galactic winds affect the transition 
from metal-free to normal star formation. Because the transition is actually spatially inhomogeneous and 
temporally extended, this mechanism cannot be responsible for double reionization given plausible parameters 
for the winds. We next consider photoheating, which causes the cosmological Jeans mass to increase in ionized 
regions and hence suppresses galaxy formation there. In this case, double reionization requires that small halos 
form stars efficiently, that the suppression from photoheating is strong relative to current expectations, and 
that ionizing photons are preferentially produced outside of previously ionized regions. Finally, we consider 
H2 photodissociation, in which the buildup of a soft ultraviolet background suppresses star formation in small 
halos. This can in principle cause the ionized fraction to temporarily decrease, but only during the earliest 
stages of reionization. Finally, we briefly consider the effects of some of these feedback mechanisms on the 
topology of reionization. 

Subject headings: cosmology: theory — galaxies: evolution — intergalactic medium 



1. INTRODUCTION 

The epoch of reionization marks the time when collapsed 
objects began to influence the diffuse intergalactic medium 
(IGM) and first rendered it transparent to ultraviolet photons. 
It is an important milestone signaling the end of the "cos- 
mological dark ages" and the emergence of the first lumi- 
nous sources in the universe. A wide variety of observational 
probes have recently been used to constrain not only the reion- 
ization process itself but also the properties of the sources 
driving it. 

The most powerful observational constraints co me from 
Lyg absorption spectra of high reds hift quasars ( Beck er et alJ 
I200U iFan et alJl2(X)l |mite et al." 2003). At least one z > 
6 quasar shows a complete Gunn & Peterson ( 1965) trough 
as well as a rapidly evolving neutral fraction indicative of 
the final stages of reionization (though see Songaila 2004 
for a different interpretation). This conclusion is strengthed 
by analyse s of the proximity effects around these quasars 
JWviflie & Loeb 2004a; Mesinaer & Haiman 2004). The sec- 
ond clue is from measurements of the large scale polariza- 
tion anisotropics of the the cosmic microwave background 
(CMB), which imply a high optical depth to elec tron scatter- 
ing and require reionization to begin at z > 14 ( Ko gut et alJ 
|2003; Spergel et al. 2003|, albeit with large error bars. A 
third clue comes from the temperature of the Lya forest at 
z ~ 3, which is relatively high and indicates that the neutral 
fraction cha nged substantially at z < 10 dTheuns et alJi2002t 
iHui & H aiman 2003). However, this constraint is difficult to 
interpre t because helium reionization also heats the IGM (e.g., 
LSokasian etalj2002ft. 
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Taken together, these constraints imply a complex and ex- 
tended reionization history, ruling out the simplest models in 
which the emissivity is a function purely of halo mass and 
independent of cosmic time (e.g., Barkana & Loeb 2001 and 
references therein). Within the context of the cold dark mat- 
ter model, measurements of the matter power spectrum imply 
that structure formation occurred extremely rapidly at high- 
redshifts: that is, the fraction of baryons incorporated into 
star-forming halos at z = 6 is over three orders of magnitude 
larger than the fraction at z = 20. If the source properties re- 
mained constant, this exponential growth of structure would 
require reionization to be completed within a short time in- 
terval. Thus, if reionization is to begin at z = 20 but end at 
z — 6, the sources themselves must evolve significantly dur- 
ing reionization. 

Such evolution could occur naturally through stellar 
feedback, because the first galaxies form out of metal- 
free gas and have exceedingly shallow gravitational po- 
tential wells. A number of models have been devel- 
oped to reconcile the disp arate i nference s about reionization 
(Wvithe&Loeb 2003a b: J c^ l2003albt iHa iman & Holdej 
12003; Sokasian et al. 2004'; 'Fukusita & Kawasaki 2001 
ISpmervifle & Livio 2003; Onken & Miralda-Escude 2004), 
all incorporating some sort of feedback mechanism to de- 
crease the ionizing efficiency of the sources. The resulting 
ionization histories display a wide range of features and can 
extend over long redshift intervals. One of the more intriguing 
possibilities is so-called double reionization. We will use this 
term to designate histories in which the global ionized frac- 
tion icj decreases with cosmic time over some interval. Such 
models provide the clearest signatures of feedback and should 
be the easiest to observe with future 21 cm tomography mea- 
surements (Furlanetto et al. 2004c). They can be contrasted 
with "stalling" models in which .f, increases monotonically 
but the feedback mechanism still shapes the overall evolution. 

In this paper, we will critically examine the plausibility of 
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the physical processes that may lead to double reionization. 
Existing models attribute double reionization to three differ- 
ent feedback mechanisms. The first is the transition from 
metal-free Population III (hereafter Pop III) star formation to 
"normal" Population II (hereafter Pop II) star formation as ob- 
served in the local universe. If Pop III stars are massive, they 
can produce approximately an order of magnitud e more ion- 
izing photons per baryon than do normal stars ( B romm et alJ 
I^Olb). The first stars therefore efficiently ionized the IGM; 
however, as Pop III stars died and exploded, they expelled 
metals and enriched the IGM. Once star-forming regions 
reached a typical metallicity of Z 10~^'^ Z©, the excess 
cooling provided by the metals could reduce the Jeans mass 
and switch the star-formation mode from P op III to Pop II 
jBromm et al. 2001a: Bromm & Loeb 2003). thus lowering 
the ionizing efficiency. Under some circumstances. Pop II 
stars could no longer counteract recombinations and so x, 
would decrease. Unfortunately, existing models have treated 
this tra nsition crudely. Both Cen (2003a) and Wvithe & Loeb 
(|2003a) assigned a single, universal redshift at which star for- 
mation switches between the two modes; the redshift was 
taken to be the time when the mean metallicity of the universe 
passed the above threshold. Both found that double reioniza- 
tion occurred for a range of input parameters. However, dif- 
ferent regions of the universe are expected to reach the same 
evolutionary stage of structure formation at different cosmic 
times due to modulation by inhomogeneities on large scales; 
this cosmic variance is particularly large at high redshifts 
(iBarkana & Loeb 2004). Haiman & Holder (2003) pointed 
out that if the Pop Ill/Pop II transition is spread over a simi- 
larly long time interval, double reionization no longer occurs. 
However, they did not calculate the plausible duration of the 
transition, and they conflated the metallicity transition with a 
change in the mass threshold of galaxy halos for star forma- 
tion. In |2j we will explicitly compute when and how rapidly 
the Pop Ill/Pop II transition occurs using a physical model for 
enrichment by galactic winds. Like Scannapieco et al. (2003), 
we find that the transition must occur over an extended red- 
shift interval. In such a scenario, double reionization requires 
a number of unlikely assumptions. 

A second feedback mechanism is photoheating. Reion- 
ization raises the IGM temperature from T < 100 K to T > 
10'* K, increasing the ambient pressure and hence the cos- 
mological Jeans mass llReesl [l986: Efstathiou KTool . As a 
result, low-mass halos can no longer collapse in ionized re- 
gions and the star formation rate may decrease sharply. The 
degree of suppression is not clear. Early work suggested 
that this mechanism prevents halos with circular velocities 
Vr 5, 30-50 km s~' from forming (Thoul & Weinbers 199S 
iKitavama & Ikeuchill2()00 ). However, Diikst ra et aL (>2004al) 
showed that the suppression is considerably weaker near the 
time of reionization, because in that case many halos have al- 
ready begun to collapse and their high densities shield them 
from the ionizing background. At the same time, exist- 
ing halos with virial temperatures Tyir < 10"* K will photo- 
evaporate as they absorb ionizing photons (Barkana & LoeQ 
[1^99; Shapiro et al. 2004). The net effect is that, as reion- 
ization proceeds, the mass threshold for galaxy formation in- 
creases by some (uncertain) amount and reionization slows 
down. Again, in such a situation it is possible for the ionizing 
sources to "overshoot" and for recombinations to dominate 
for a time. This mechanism has been included in most exist- 
ing models and rarely been found to produce double reioniza- 
tion. Nevertheless, there has been no systematic study of the 



requirements for such a phase. In ^ we show that photo- 
heating can in principle cause double reionization, but only in 
exceptionally optimistic circumstances. 

A third and final feedback mechanism is the photo- 
dissociation of H2. Rotational transitions of H2 provide 
a cooling channel that operat es in halos w ith Tyir > 200 K 
(jHaiman et al. 1996; Tegmar k et alJ fT997h . well below the 
threshold at which atomic cooling becomes efficient. The 
first halos to form stars in any hierarchical model have small 
masses and therefore rely on H2 for their cooling (Abel et alj 
2002; Bromm et al. 2002). However, H2 is fragile and is eas- 
ily dissociated by soft UV photons in th e Lyman- Werner band 
(11.26-13.6 eV) (iHaimanet al.lfT997l) . Once the first stars 
build up a sufficient UV background, this cooling channel ter- 
minates and the minimum halo mass to form stars increases; 
operationally this mechanism is similar to photoheating al- 
though it operates at a lower halo mass scale. We will con- 
sider whether this scenario can lead to double reionization in 

Of course, by suppressing the ionizing efficiency in biased 
regions these feedback mechanisms affect not only the global 
rei onization history but also its topology. As recently stressed 
by 'Furianetto et al.' ('2004b, hereafter FZH04), the topol- 
ogy of reionization is observable through 21 cm tomography 
(see also Furianetto et al. 2004c) and through Lya absorption 
(Furianetto et al. 2004a; Wvithe & Loeb 2004b). Although a 
much more difficult problem to approach analytically than the 
global evolution, it is also worth considering what signatures 
one might expect to see in the bubble size distribution from 
these feedback mechanisms. We briefly consider this ques- 
tion in Island then discuss other consequences of our results 
in^ 

Throughout our discussion we assume a cosmology with 
n„, = 0.3, Q-A = 0.7, Qt = 0.046, Ho = 100/1 km s"' Mpc"' 
(with /; = 0.7), « = 1, and cg = 0.9, consistent with the most 
recent measurements ( Spergel et al. 2003.) . We will use co- 
moving units unless otherwise specified. 

2. A METALLICITY TRANSITION 

As noted above, existing models attributing double reion- 
ization to the evolving cosmic metallicity rely on an instanta- 
neous transition between the two modes. In reality, there are 
two mechanisms that will slow it down: the formation of stars 
in pre-enriched halos (even early on) and the finite velocity of 
winds. Here we will examine these two processes in turn and 
apply them to the reionization history. 

2.1. The Formation Rate of New Halos 

We divide newly-collapsed gas (i.e., baryons that have 
just been incorporated into galaxies able to form stars) 
into two components. First, some fraction of the gas ac- 
cretes onto existing halos that have already formed stars. 
Provided that an existing halo had already enriched it- 
self beyond the threshold metallicity Z, to form Pop II 
stars, this component of newly collapsing gas cannot make 
Pop III stars. The metallicity thre shold is Z, ^ 10^^^ Zq 
dBromm et al. 2001a: IBromm & Loebll2003h . which requires 
A > 10-^^(0.5//z)(Z,/10-3-5 Zq), where is the star for- 
mation efficiency and fz is th e fraction of the stellar mass 
converted into heavy elements. 'Heger & Wooslev'('2002) find 
fz ^ 0.5 for stars which undergo pair-instability supernovae, 
although not all massive stars actually explode so this is only 
an upper limit to the true value. The fraction is somewhat 
smaller (but still > 0.1) for normal initial mass functions 
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(IMFs). The metallicity threshold for Pop II star formation is 
sufficiently small that we expect any halo that has previously 
form stars to make only Pop II stars after the initial starburst, 
unless the metals remain highly clumped in the halo. Because 
clumping would only slow down the transition between the 
two star formation modes, we will conservatively ignore it. 

To form Pop III stars, gas must therefore collapse into a 
"new" or "fresh" halo that is forming stars for the first time. 
We take this criterion to be a minimum galaxy mass m^in de- 
termined by the physics of radiative cooling. We will typi- 
cally assume Wmin corresponds to Tvir = 10** K, the threshold 
temperatu re above which atomic hydrogen cooling is effec- 
tive (e.g., Barkana&Loebll2001 l). Thus, to find the maximal 
rate at which Pop III stars can form, we wish to know the rate 
at which halos of mass m w m^in form. 

Unfortunately, as described in the Appendix, computing 
this rate is beyond the capability of existing semi-analytic 
models. We there fore make a simple estimate ('Sasak?'1994l; 
IVerde et alJl2()()l . Given the[P ress & Schechter ( 1974) mass 
function, we can compute the total rate of change of the halo 
mass function: 
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where v — 5c{z) / (j{m), 5c{z) is the critical density for col- 
lapse, and CT^ (m) is the variance of the density field smoothed 
on the mass scale m. This total rate of change is the difference 
between the creation and destruction rates of halos. The two 
terms in the derivative can be conveniently (if only approxi- 
mately) identified as these two rates. For rare, massive halos 
(with 3> 1), the creation term should dominate because there 
are only a small number of larger objects into which the halo 
can be absorbed. The creation rate should therefore increase 
with increasing v. Halos with 1, on the other hand, are al- 
ready common. As the nonlinear mass scale rises, these small 
halos must accrete onto more massive objects. Thus we ex- 
pect the destruction term to increase with decreasing v. The 
two terms in equation have the correct limiting behavior, 
although of course the split does not capture all of the physics, 
particularly when v ^ \ . The rate at which genuinely new ha- 
los form is approximately 
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where p is the mean cosmic density. 

The top panel of Figure [2 shows the evolution of d/„e 
relative to the total rate of collapse d/coii/dz, where 

5c{z) 



erfc 



\/2cr„ 
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v/dz 



(3) 



is the fraction of matter bound to halos with m > m^in and 
cTmin = o-(mn,in). We require that d/ne„/dz < d/^oii/dz; the ap- 
proximations in equation (|2j violate this requirement at suf- 
ficiently high redshifts. The solid line associates nimin with 
a virial temperature Tyii- = lO'* K. The dotted and dashed 
lines assume values of mmin an order of magnitude larger and 
smaller than this value, respectively. The figure reveals two 
important points. First, d/new/dz evolves slowly with redshift. 
It does not approach zero until z ~ 6, by which point reion- 
ization is complete. This is because m„^in is either larger than 
or close to the nonlinear mass scale at these early redshifts, so 
pristine galaxies are still forming. In fact, the magnitude of 
d/new/dz actually increases with cosmic time; its relative im- 
portance decreases only because the number of halos larger 
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Fig. 1 . — Evolution of the fraction of gas collapsing into new halos. Top: 
Approximate model based on eq. Bottom: Extended Press-Schechter 
model (see Appendix). In each panel, the solid lines have our default m,„i„ 
(corresponding to T^-„ = lO"* K). The dashed and dotted curves correspond to 
decreasing or increasing m,„i„ by an order of magnitude, respectively. 



than the threshold increases even more rapidly (as they corre- 
spond to a higher value of v), and each of these objects con- 
tinues to accrete mass. Second, the fraction of mass entering 
new objects increases with increasing m^^^ because such ob- 
jects become rarer and it takes the nonlinear mass scale longer 
to reach a larger threshold. On the other hand, d/new/dz de- 
creases somewhat more rapidly in this case. 

Because equation Q is only approximate, we describe an 
analogous calculation with the extended P ress-S chechter for- 
maHsm (Bondetal. 1991; Lacev & Col3 ll993h in the Ap- 
pendix. We show the results in the bottom panel of Figured] 
in this calculation the tran sition b etween the two modes is 
even slower (see also Sca nnapieco e t al. 2003). Thus a slow 
transition seems inevitable, and we will clearly require winds 
to introduce sharp features in the reionization history. 

2.2. Wind Model 

Galactic winds are the m ost likely agent for spreading met- 
als through the IGM (e.g., 'Aauirre et al."200r, 'Madau et alj 
BoOl; Scannapieco et al. 2002; Furlanetto & Loeb 2003) and 
for accelerating the Pop Ill/Pop II transition. In this section 
we will outline a model describing them. In order to keep 
our methods as straightforward as possible, we will examine 
a simple model for the wind expansion that captures the ex- 
pected scaling laws and calibrate it against more advanced 
treatments. 

The key input is the underlying energy budget of the 
winds. The available energy from supernovae is Wsn = 
fi,EsNmh/oJsN, where is the star formation efficiency 
(which may depend on halo mass), Esn — lO^'/isi ergs is the 
energy per supernova, ujsn is the mass in stars per supernova 
event, and m/, is the baryonic halo mass. We fix Esn and ujsn 
here for simplicity. For Pop II stars, we tak e £51 = 1 and 
LOsN = 100 M0, in accord with a lScalol(ll998h IMF for stellar 
masses in the range 0.1-100 M0. For Pop III stars, we take 
£51 = 10 and u!sN = 100 M0, assuming that most of the stel- 
lar mass undergoes pair-instability supern ovae and maximiz- 
ing enrichment from the first generation jHeeer & WooslevI 
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l2(K)l . 

The simplest possible wind model assumes an energy- 
conserving point explosion in a constant density medium. 
This well-known Sedov (1959) solution predicts that R cx 
{WsNt^/pY^^, where R is the radius, p is the density, and t is 
the elapsed time. Although this solution neglects a number of 
important processes (see ^2.2. 1> . iFurlanetto & Loebl (120031) 
showed that the scaling remains approximately correct, al- 
though the normalization does change significantly. Then we 
can express 77, the ratio of total mass enriched by the wind to 
the mass of each galaxy, as 

where K„ is a normalization constant that accounts for the 
many factors we have neglected (note that it could also in- 
clude variations in ujsn)-^ We have assumed that f w 1 /3H{z) 
(i.e., each wind has propagated for half of the age of the uni - 
verse). Comparing to the results of Furlanetto & LoebH2003l) . 
this solution yields radii 2-3 times too large, depending on 
the mass and redshift, suggesting that typical wind models 

have K„, ^ ^-^f^/^, where /,v accounts for radiative losses. 
Of all possible cooling mechanisms, we include only Comp- 
ton cooHng by setting /h, min(l,fcomp /?//), where fcomp is 
the Compton cooling time and w i/^' (z). We ignore other 
possibilities (most importantly those in the dense interstellar 
medium of the host) in order to maximize enrichment. 

Note that the two stellar populations have different efficien- 
cies 77. In evaluating equation for a particular galaxy, we 
assume that the fraction of Pop III stars is /new//coib where 
/new is the total mass fraction in Pop III stars integrated over 
time. This does slightly underestimate the effects of Pop III 
stars because smaller halos have more Pop III stars and are 
more efficient polluters. However, a more detailed treatment 
would require fully self-consistent merger rates, which are 
currently unavailable. 

2.2.1. Caveats 

Our simple model neglects a number of important effects, 
which we list here for completeness. We emphasize that 
[Furlanetto & Loeb (2003) included all of these and showed 
that they do not affect the scaling. The first is radiative 
cooling, which we have already discussed. Second, the Se- 
dov solution neglects the gravitational potential of the host 
halo. At large masses, gravity tends to flatten the mass de- 
pendence (Furlanetto & Loeb 2003); however, such halos are 
rare at these high redshifts. A third difficulty is in the choice 
of elapsed time, because halos have a range of star forma- 
tion histories. Crudely, we would expect more massive ha- 
los to be older, which steepens the mass dependence. On 
the other hand, the Sedov scaling also neglects the evolving 
background density (and the halo density profile). Older halos 
spend more time expanding into dense environments, which 
flattens the mass dependence. Finally, we ignored the fact that 
the wind travels through an expanding medium, so that it need 
n ot accelerate swept u p material from rest. 

IFurlanetto & Loebl 0200 1) noted that we can set an upper 
limit to the wind size by balancing the energy input with the 
energy required to accelerate the enclosed baryons to the Hub- 
ble velocity at the outer edge of the wind. This procedure 

' In the following, we actually use 1 + r]{m) to include self-enrichment of 
the galaxy. 



yields rj^.^^ w l.li] with exactly the same scalings as in equa- 
tion 0. Thus the above estimate is reasonable; it is signifi- 
cantly smaller than the maximal limit because the winds have 
not completed their expansion. 

2.3. The Enrichment History 

We now wish to compute the probability that a collapsing 
halo (with m w m^in) forms in a region already enriched by 
galactic winds. As a first step, given the enrichment efficiency 
ri{m), we can estimate the fraction of space with metals: 

e;(z) = ^ dm (^-^^ r]{m)n{m). (5) 

This equation would be accurate if the galactic winds did not 
overlap or if their volumes were additive (as is the case with 
H II regions, which occupy a net volume dictated by the total 
number of ionizing photons). Because winds expand at much 
less than the speed of light, the latter is not a good approxi- 
mation. If the host galaxies were randomly distributed, then 
the true filling factor would be = 1 — exp{—Q'^). 

However, rather than the simple volume fraction of en- 
riched material, we actually wish to compute the probability 
that a new halo forms in an enriched region . Collapsing ha- 
los are biased and therefore lie near existing halos (and their 
winds). We will now estimate the importance of this effect. 
We first define R„ to be the average wind size (in comoving 
units). The excess probability that two galaxies sit near each 
other is parameterized by the correlation function £^gg, which 
to linear order can be written as £1,^ = bib2^s5, where bi and 
b2 are the halo biases (Mo & White 1996) and ^ss is the dark 
matter correlation function. For simplicity, we will use the 
lin ear power spectrum for the latter, with the transfer function 
of Eisenstein & Hu (1998). 

The probability that a newly-forming halo lies within R„ of 
an existing galaxy is then approximately 

Qe = Q'Ai+b,,^b,4s5iR.)l (6) 

where b^ is the mean bias of a wind, and bmw is the bias of 
a newly-formed halo.** The mean bias of the enriched regions 
is: 

J dmmri(m)b{m)n{m) 
bw = n-: 1—^ , (7) 

where the extra factor of m enters because we compute the 
bias of the enriched volumes rather than the bias of ha- 
los. The bias of newly-formed galaxies is b^^^ w ^('Wmin) 
because all such objects have mass close to the low-mass 
threshold for star formation. Finally, we assume randomly 
distributed wind hosts by setting the enrichment probability 
Pe = 1 -exp(-g^). 
Figure|2]shows the resulting enrichment histories. We take 
= /™ = 0.1; the different curves take different wind ef- 
ficiency parameters K^.. The thick curves in the top panel 
show /^pristine = ^ — Pe, the probability that a new halo forms 
Pop III stars. The time at which Pop III stars disappear ob- 
viously depends strongly on the wind efficiency parameter. 
For realistic choices {K^ ^1), pristine halos become rare 

1 /3 

at z < 10. The thin curves (shown only in the Kj =1/3 
and 1 cases) illustrate the behavior if we neglect cluster- 
ing: because galaxies are so highly biased at these large red- 
shifts, clustering makes enrichment of galaxies much easier. 

Taking the average i^gg over the wind volume does not significantly affect 
our results. 
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Fig. 2. — Top: Probability tliat a new halo forms out of pristine gas. The 
thick curves include clustering, while the thin curves neglect it. Bottom: The 
total fraction of collapsing gas forming Pop III stars; all curves include clus- 
tering. We take /|' = /|" = 0.1. The solid, long-dashed, short-dashed, and 
dotted curves assume kI/^ = 1/3, 1/2, 2/3, and 1, respectively. 



It completes when only about half of the total volume has 
been injected with metals. The bottom panel shows the the 
total fraction of collapsing gas able to form Pop III stars, 
^pristine X (d/new /dz) / (d/coii /dz) . Although Q'^ evolves fairly 
rapidly, including the fraction of gas that accretes onto exist- 
ing halos significantly moderates the rate of decline. In this 
model, the evolution of d/new/dz dominates the Pop Ill/Pop II 
transition until relatively late, regardless of K„. 

2.3.1. Caveats 

We now outline and justify the simplifications that led to 
equation (|6}. The most serious shortcoming is our treatment 
of wind overlap. We have made no attempt to follow the 
dynamics of colliding winds. However, since the wind vol- 
ume scales with the deposited energy more slowly than lin- 
early, combining their energies through overlap will only re- 
duce the total enriched volume relative to a case without over- 
lap. Thus the net effect of overlap is to slow the transition 
from Pop III t o Pop II and therefore strengthen our conclu- 
sions in ^2.41 Second, we have assumed that galaxy for- 
mation in the wind-enriched regions can be well-described 
by dark matter dynamics. This may not be a good assump- 
tion if the winds disturb a large fraction of the IGM by, for 
example, sweeping it into a shell, stripping nearby halos as 
they collapse, or triggering sta r formation behind the termi- 
nation shock (I Scannapieco et al. .2001; Cen ..2003a) . If the 
feedback is positive, double reionization would become even 
more difficult, but if the winds suppress galaxy formation the 
prospects for double reionization would be enhanced. (That 
case would be qualitatively similar to photoheating, except 
that the feedback would operate much more slowly.) 

There is also a set of issues related to our simplified global 
treatment of enrichment. Equation (|6j only includes linear 
halo bias, and we only use the linear ^55. The latter is proba- 
bly not a bad approximation, because small-scale nonlineari- 
ties are primarily due to the density profiles of i ndividual ha- 
los (e.g., the halo model: ICoorav & Shethl2002h . In our case. 



all the baryons inside each halo likely collect into a single 
galaxy , in which case the one-h alo nonlinearities are not rel- 
evant. Scanna pieco et alJ (|2003) constructed a more detailed 
enrichment model that include d nonlinear bias. They used the 
two-point halo mass function ( Scannapieco & Barkana 2003) 
to compute the probability that a halo lies within the winds 
of its neighbors. Our extended enrichment histories are qual- 
itatively consistent with theirs, so the nonlinear corrections 
included by this formalism do not seem to affect our conclu- 
sions. Another crucial simplification is that we have taken the 
wind sources to be randomly distributed. In reality, the bias of 
the wind hosts will increase the amount of overlap, decrease 
Qe, and extend the transition; in this sense our conclusions 
about double reionization are conservative. All of these is- 
sues are best addressed with numerical simulations, which we 
forego in order to keep our models as simple as possible. 

2.4. The Reionization History 

We are now in a position to determine how enrichment, and 
the accompanying transition from Pop III to Pop II star for- 
mation, affects reionization. The rates at which Pop III and 
Pop II stars form are 



d/m _ , > d/nev 

-^-PpristinelZj-^ 



and 



d/n 
dz 



d/cc 



dz 



d/m 
dz ■ 



(8) 



(9) 



To translate these into production rates of ionizing photons, 
we define Q as the number of ionizing photons produced per 
collapsed baryon for Pop / (where / =11 or III), e.g. Qi = 
^Hefifesc^'-yh' ^i*-^ /esc the fraction of ionizing photons able 
to escape the host halo, Nji, the number of ionizing photons 
produced per baryon incorporated into stars, and Ahc a cor- 
rection for helium. We will fix A^;^'^ = 4000 (appropriate for a 
Salpeter I MF) and N^l = 30,000 (appropriate for very mas- 
sive stars: iBromm et al. 2001b). Note that affects both 77 
and C, while /esc and K^. act as separate efficiency parameters 
for the two mechanisms. At least crudely, these are the source 
parameters that we can tune in order to test the plausibility of 
double reionization. 

The number of ionizing photons emitted per baryon per unit 
redshift is then 



. ^ , d/i d/u 
e(z) = Cn-5- + Ciii-5— : 
dz dz 



(10) 



and t he global ionized fraction x, evolves a ccording to 
( Sh apiro &Girouxl 19871 iBarkana & Loebl2001h 



dx, 
dF 



: e(z) ~XiaBneC{l +z) 



(11) 



where is the case-B recombination coefficient, n^ is the 
comoving electron density in a fully ionized medium, and C 
is the density clumping factor within ionized regions. The 
first term is the rate at which stars produce ionizing photons, 
while the second is the rate at which protons and electrons 
recombine. We treat C as a constant; a better model would 
include t he evolving density d istribution of ionized gas in 
the IGM. iMiralda-Escude etal.. (2000) have developed such 
a model calibrated to numerical simulations. They argue that 
low-density regions, where the recombination rate is small, 
will be ionized first. As the ionized volume increases, the ion- 
izing background can keep progressively denser regions ion- 
ized. They showed that while C increases throughout reion- 
ization, it does so only modestly until x, > 0.9, because nearly 
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Fig. 3. — Reionization histories, (a): Normalized emissivity in our fiducial 
model (with /J' = = 0. 1 and /''^ = /™ = 0.05). The solid, long-dashed, 

1 /3 

short-dashed, and dotted curves assume A",/ = 1/3, 1/2, 2/3, and 1, re- 
spectively, (b): Same as (a), except /'JJ, = 0.4. (c): Ionized fraction in the 
fiducial model, assuming C = 1. (d): Ionized fraction for the models in (h). 
All curves assume C = 3 except the dot-dashed line, which has C = 1 but 
is otherwise the same as the dotted curve. The corresponding enrichment 
histoiies are shown in Figurel2l 



all of the volume is near the mean density at the high red- 
shifts of interest. Since we are primarily interested in the 
middle phases of reionization, the assumption of a constant 
C is reasonable for most of reionization unless there is con- 
siderably more small-scale structure than the si mulations o f 
iMiralda-Escude et al. (2000) found (e.g., Haiman et all2001l) . 
However, it will underestimate the recombination rate at the 
end of reionization. (Note also that with a constant C we must 
explicitly forces, < 1.) 

Figure |3plots the reionization history in two sets of mod- 
els. The top panels show the evolution of e{z) and the bot- 
tom panels show We vary the wind efficiency from 

kI^^ = 1 to = 1/3 in each panel. In all cases, we as- 
sume = /™ = 0.1 as in Figure |2j because these quanti- 
ties affect both enrichment and reionization they have little 
effect on our conclusions. Panels (a) and (c) show a base- 
line case in which fH^. — /^'J, ~ 0.05. The difference in emis- 
sivity between the two generations is then determined purely 
through the stellar physics encapsulated by A^.^^. In panels (b) 
and (d), we set flf^. — 0.4 in order to exaggerate the differ- 
ence between the generations (although there is some moti- 
vation for a high escape fraction from the first generation of 
sources; ^Whalen et al. 2004). In panel (c), we assume that 
C = 1, as is appropriate for the bulk of the IGM according to 
IMiralda-Escude et ab (12000)) . In panel (d), however, we as- 
sume C = 3. 

According to panel (a), the difference in A^-^^ between Pop II 
and Pop III is clearly insufficient for double reionization: the 
slow evolution in d/new/ dz smooths out the expected discon- 
tinuity. Not surprisingly, is therefore monotonic regardless 
of the wind efficiency. Interestingly, this is true even though 
we have "tuned" the parameters so as to force complete en- 
richment and reionization to occur at approximately the same 
epoch. Thus, to introduce features into x, we must exaggerate 



the contrast between the generations, as in panels (b) and (d). 
By manipulating the escape fractions, we cause a clear drop in 
the emissivity during the Pop Ill/Pop II transition. However, 
even in this case we still only force double reionization if the 
winds are powerful and clumping is strong. We remind the 
reader that even with f„, — 1, the wind enrichment models of 
l^rlanetto & Loeb (2003) had K,, < 1 /8, so we regard larger 
values as implausible. Double reionization is difficult because 
maintaining full ionization requires only that the production 
rate of ionizing photons exceeds the rate at which hydrogen 
atoms recombine (i.e., the right-hand side of equation dl 1> 
must be positive). Thus Xj will only decrease with time if 

e(z)<0.1Cx-- f^j . (12) 

The marked decrease for the solid and long-dashed curves is 
large in relative terms but still leaves the absolute emissiv- 
ity strong enough to cancel recombinations. This illustrates 
an obvious feature of these models: if reionization completes 
well before enrichment, the emissivity will have risen suffi- 
ciently between these two epochs so that, even after a sharp 
fall in e, the universe remains completely ionized. If we ac- 
celerate enrichment through more powerful winds (the short- 
dashed and dotted curves), e does fall below the threshold of 
equation il2\ . and a brief epoch of recombination can ensue. 
However, this is only true if C is larger than unity. The dot- 
dashed curve shows x, for the same model as the dotted curve 
(i.e., maximal winds) but with C = 1. In this case, x,- flattens 
but does not decline significantly. 

We have explored a range of parameter choices and found 
that double reionization can in principle be achieved, but only 
in rather extreme models. In order for the emissivity evolution 
to have a strong feature, we require Cm ^ Cu (and in partic- 
ular much larger than the expected difference in Nyt>). Sec- 
ond, enrichment must be timed to occur during or soon after 
reionization, so that the emissivity has not risen too far before 
enrichment is complete. Third, the magnitude of the emis- 
sivity after the Pop Ill/Pop II transition must satisfy equation 
il2\ . This is easiest to achieve by boosting the IGM clump- 
ing factor, because that accelerates recombinations and gives 
the emissivity less time to recover. We find that double reion- 
ization only occurs in parts of the IGM with C > 3. Because 
such regions fill only a small fraction of the volume in stan- 
dard models (e.g., Miralda-Escude et aL.2000.) . a recombina- 
tion phase is strongly disfavored. 

2.4.1. Quasars and a Variable f-^ 

There is some evidence that cx cr^ oc rn^/^ ( 1 + in lo- 
cal galaxies with velocity dispersions a < <Jr = 124 km s~^ 
(Dekel&Woo 2003; Kauffmann et al. 120031) : above this 
threshold the star formation efficiency remains constant. How 
would such behavior in high-redshift galaxies affect our con- 
clusions? We let = f*{pc) and assume f^, cx m^/^ for all 
masses. (Massive halos with a > Oc are so rare at high red- 
shifts that they have a negligible effect on our results.) Note 
that the Pop III star formation efficiency is a constant because 
such stars form only in halos with m ^ mmin. We can repeat 
our reionization calculations using this prescription; the only 
adjustment is that we must modify equation dlOt to account 
for the increasing as galaxies grow (e.g., Wvithe & Loe^ 
'2003a). For a fixed energy input, the fi, cx m^/^ models tend 
to have somewhat slower enrichment because r\ is suppressed 
in the (abundant) low-mass halos in which winds are most 
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efficient (eq. 0). It is also somewhat harder to achieve a 
sharp drop in the emissivity because the typical galaxy mass 
increases with cosmic time (and hence so does the mean 
However, the law has no real effect on our conclusions re- 
garding double reionization: our direct calculations for the 
case of cx imply that the same fairly extreme conditions 
outlined above must be fulfilled. Thus, for the sake of brevity, 
we have chosen not to show any explicit results for these mod- 
els. 

We note here one interesting result of our models: if cx 
m^/^, early reionization at z 3> 6 requires the Pop III star 
formation efficiency to be much larger than the Pop II effi- 
ciency in halos of a comparable mass. The reason is simply 
that Pop III stars must form in pristine halos with m = m„iin, 
where fl^ « 0.01/JJ.. The number of Pop III stars would there- 
fore also be extremely small, unless fi"(mn,in) ^ fl^{mj„i„). 
Previous studies, such as|iWvithe & Loeb (2003a), had found 
such early reionization to be possible because they allowed 
more massive halos (where the star formation efficiency is 
large) to form Pop III stars even though such galaxies would 
have been self-enriched. Thus, if observations continue to fa- 
vor early reionization, it appears the star formation efficiency 
must scale differently with mass than in nearby galaxies, at 
least for one of the two star formation modes. 

Similar models could also describe the i nput of black holes, 
which may generate powerful winds (see |Furlanetto & Loebl 
120011 and references therein). Wvithe & Loeb (200^ find 
a good fit to the quasar luminosity function over the range 
z = 2-6by assuming that Mbh ~ 10^'' M0(fT/54km s"')^. We 
can estimate how quasars affect double reionization by as- 
suming that this relation extends to the higher redshifts and 
lower masses of interest. We will also assume that a quasar 
radiates a fraction eg of its rest mass and that a fraction <7vv 
powers the wind. So long as ^constant, quasars make only 
a small difference to the results if eg ~ 0. 1 because black 
holes can be ignored in the small halos responsible for most 
of the enrichment and ionizing photons. On the other hand, if 

cx m^/^, quasars can become significant in small halos. For 
example, the ratio of wind input energies is 



Wq ~ ■ ^0.1 q,, tQ " 



10 



5/2 



(13) 



Quasars can be inserted directly into our cx m^/^ models 
without modification; clearly they will not alter the likeli- 
hood of double reionization, although they could alleviate 
the problem mentioned in the previous paragraph (but see 
iDiikstra et al..2004b) . 

3. PHOTOIONIZATION HEATING 

The second type of feedback is radiative, for which the 
most important example is photoionization heating. As de- 
scribed in ^ this raises the Jeans mass in ionized regions by 
an uncertain amount and suppresses accretion onto existing 
small halos. Our primary interest is in whether this feedback 
mechanism can cause double reionization; we will therefore 
take a phenomenological approach and leave the minimum 
galaxy mass in heated regions as a free parameter. Because 
the physics depends on the depth of the potential well, we 
will phrase the threshold in terms of a (redshift-independent) 
virial temperature T/,. We will also let = rvii-(w^min); the ^ 
and c subscripts refer to hot and cold gas, respectively. We 
will take T/, = 2.5 x 10^ K as a fiducial value, corresponding 
to w 50 km s^'. Once a region is ionized, it will begin to 



cool adiabatically (and radiatively if it recombines). This will 
soften the increase in the time-averaged Je ans mass (or "filter 
mass ': lGnedin & Huai99a) . ICer] ll20033 ) has included these 
effects more precisely through a "phase-space" description. 
By ignoring this possibility, we will maximize the possibili- 
ties for double reionization. 

A number of authors have already included pho- 
toionization heating in models of reionization (e.g^ 
Wvithe & Loeb 2003a; Cen 2003a; Somerville & Livi3 
2003; Onken & Miralda- Escuda 120041) . In this case the 
emissivity has two terms: 



e(z) = Xh 



dz 



dz 



(14) 



where x/, is the filling factor of regions that have been heated, 
fh = /coii(> 7},), and /,. = /coii(> 7^). We take the deriva- 
tive of the mass-averaged ionizing rate for each component, 
(C//,) and (C/f), rather than of the simple collapse fractions 
to include the possibility that the star formation efficiency is a 
function of halo mass. (The functional dependence of qual- 
itatively affects the results for photoionization heating, unlike 
in the metal enrichment scenario.) We may then insert this 
emissivity into equation (II 1> and solve forx, (z). Note that we 
assume both kinds of sources to have the same C,. 

However, to do so we must specify how Xft(z) evolves. The 
simplest assumption, and one often used in the literature, is to 
set Xi = Xh'- the fraction of gas that has been heated is the same 
as the ionized fraction. If this prescription is accurate, dou- 
ble reionization from radiative feedback is impossible even in 
principle, because x,- regulates the feedback and there is no 
possibility of "overshoot." However, we must actually have 
Xh > Xi because some photons are lost to recombinations. For 
example, suppose T/, ^ Tc, so that the emissivity in hot re- 
gions is negligible. The first halos (with Tyir — T^) collapse 
and ionize their surroundings; structure formation effectively 
stops in these ionized bubbles. However, small halos can con- 
tinue to form in the neutral gas outside these bubbles. Thus 
ionizing photons will preferentially appear in cold regions 
while recombinations will only occur in hot regions. The fill- 
ing factor of hot gas must therefore evolve faster than x,-. The 
extreme limit would be to calculate x/, without including any 
recombinations. 



dxh , ^ 

dF = ^(^)- 



(15) 



Of course, reality lies somewhere between these two ex- 
tremes: some fraction of the ionizing photons are produced 
in pre-ionized regions and counteract recombinations before 
ionizing new material, while recombinations also occur in the 
newly-formed ionized bubbles. The complexity of radiative 
transfer and feedback makes a self-consistent choice difficult 
and is probably best addressed with simulations. We choose 
equation (I15> in order to maximize the possibilities for dou- 
ble reionization; even in this extreme case we will find double 
reionization to be quite difficult. 

Figures |^ and c show reionization histories for a set of 
models with C = 100 and =constant (note that we only con- 
sider a single stellar population here). The solid, long-dashed, 
short-dashed, and dotted curves have Th = 10^, 10^, 2.5 x 10^, 
and 10^ K, respectively, with C—\. The dot-dashed curve has 
Th = 2.5 X 10^ K and C = 3. In this case, the emissivity is the 
same as for the short-dashed curve. The solid curve shows 
the ionization history without photoheating; it rises steeply, 
completing reionization in a relatively short time period. An 
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Fig. 4. — Same as Fig.|3] but for models in which photoionization heating 
raises the mass tlireshold for collapse. Panels (a) and (c) assume =constant 
and ^ = 100, while panels (b) and (d) assume /* oc m^^^ and = 500. The 
solid, long-dashed, short-dashed, and dotted curves haveT;, = 10^*, 10',2.5x 
10^, and 10* K, respectively, with C = 1. The dot-dashed curves have T;, = 
2.5 X 10^ K and C = 3; in this case the emissivities are the same as the short- 
dashed curves. 



increase in Th/Tc introduces a stronger break into the emis- 
sivity. As long as xi, <C 1, small halos dominate the ioniz- 
ing photon budget and the evolution is nearly independent of 
Th/Tc- Once x/, becomes large, e cx d(C//,)/dz and 7), deter- 
mines the amplitude as well as the duration of the break. Note 
that the redshift at which jc/, — 1 corresponds to the minimum 
in e. Unlike in the metal enrichment case, this must occur 
while Xi <\. Thus, any feature in the ionization history due 
to photoionization will occur before reionization is complete. 

For this reason, a sharp drop in the emissivity is more diffi- 
cult to achieve, because xi (or in our model the closely related 
quantity Xh) regulates the emissivity and causes it to flatten 
substantially before it can drop. Equation M2\ again expresses 
the criterion for double reionization: the emissivity must fall 
below the recombination rate. Not surprisingly, it therefore 
requires either extremely large suppression (7), ~ 10^ K) or 
substantial clumping. Our fiducial choice T^ = 2.5 x 10^ K 
is just on the threshold for achieving double reionization if 
C= 1. 

Figures EJ? and d show reionization histories for models 
with oc nP'/^ and C,c — 500 (evaluated at the fiducial Uc). 
In this case the behavior is obviously quite different, with x,- 
remaining monotonic in all cases. Even for the most extreme 
choice of T^jT^ e(z) exhibits only a slight turnover Extended 
reionization histories are more difficult to achieve because the 
star formation efficiency in low-mass halos is small, making 
their contribution to the global emissivity low and their sup- 
pression much less significant. The decline in emissivity is 
much less severe in this case, and double reionization requires 
truly extreme parameters. This also implies that if quasars 
are important in reionization and Mbh c>c cr^, they will tend 
to wash out double reionization because massive halos con- 
tribute a large fraction of the ionizations. 

We have not included halo bias in the feedback treatment in 
this section. This is unlikely to affect our final results for two 
reasons. First, the bias should be less important than for metal 



enrichment, because overlapping H II regions grow quite large 
even during the middle stages of reionization (FZH04), span- 
ning several comoving Mpc when x, ^ 0.5. On these scales 
is small and the ionized regions are probably close to a 
fair sample of the universe. Because of this overlap, it is also 
extremely difficult to calculate the relevant scale from first 
principles (see ^for a first look at the issues involved). The 
second reason is that photoheating guarantees that bias will 
modulate ionizations and feedback in the same way. In other 
words, if halos are overabundant in ionized regions, the emis- 
sivity will begin to decline at a smaller x,-. But both the heating 
and ionizing rates will decline together, so the net effect will 
simply be for feedback to flatten the emissivity sooner. The 
bias will not be completely degenerate if we use equation ( I15> 
but is still not likely to have a strong effect because that sce- 
nario essentially assumes that all ionizing photons come from 
widely separated (and hence unbiased) sources. In any case, a 
full treatment of the interplay between bias and feedback will 
have to await self-consistent numerical simulations. 

In summary, photoheating can in principle cause double 
reionization, but only under the following restricted condi- 
tions: (/) the suppression of structure formation is more se- 
vere than current estimates (T;, 3> 7"^) and/or the universe is 
clumpy, star formation is efficient in small halos, and 
{ Hi) nearly all ionizing photons go into ionizing new material, 
rather than counteracting recombinations in ionized regions 
(so that x/j > X,). Taken together, we regard these conditions 
as implausible, and conclude that double reionization is un- 
likely to occur through photoheating. 

4. PHOTODISSOCIATION OF H2 

For completeness, we now examine another radiative feed- 
back mechanism. As described briefly in ^ gas in the 
first collapsed halos cools through rotation al transitions of 
H2 JHaiman etalJII996t iTe gmark et aD ^7). which oper- 
ate for Tvir > 200 K. However, such molecules are frag- 
ile and a weak soft UV background of ^21 ^ I0~^ suffices 
to dissociate them, where = J2\ x erg cm ~^ s~' Hz' sr^' 
(Haimanetal. 2000; Machacek et aLl l200Ih . Crucially, any 
photon with energy above 11.26 eV can dissociate H2, so 
(unlike in the case of ionizing photons) any individual halo 
sees Lyman-Werner sources out to cosmological distances 
jHaiman et alJl997h . The UV background redward of the Ly- 
man limit builds up quickly as stars form, suppressing cooling 
via H2 and increasing the minimum halo mass for star for- 
mation from Tc ^ 200 K to the atomic cooling threshold of 
Th = lO"* K. 

Because the effects of H2 photodissociation can be phrased 
in a similar way to photoheating (except that the "heated" vol- 
ume is automatically larger than the ionized volume), we can 
again find scenarios in which x, decreases with cosmic time, 
as long as the clumping factor is substantial and T/jTc ^ 1. 
However, it is easy to show that this phase occurs early in the 
reionization process, because the photodissociation thresh- 
old must be reached well before reionization is completed 
(Haiman et al. 2000). We wish to estimate in the Lyman- 
Werner bands; neglecting the line opacity of the IGM, we 
have 

l+z 





dt 




dz' 



eLw(i^') 



1 



(16) 



where is the emissivity per unit frequency, v' ^ v{\ + 
z')/(l and the last factor accounts for the cosmologi- 
cal redshifting of the photon energy. The upper limit Zmax 
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enters because ionizing photons cannot propagate through the 
neutral IGM. The Lyman s eries introduces a " sawtooth modu- 
lation" of the background jHaiman et alJll997 ). but this does 
not affect our estimate. To connect to the ionized fraction, 
we write euwii^) =hv y. {\/v)y. xnion, where nio,, is the rate 
at which ionizing photons are produced per unit volume and 
X is the number of Lyman-Werner photons per frequency 
decade divided by the number of ionizing photons per fre- 
quency decade.^ Ciardi & Madau (2003) have evaluated a 
closely related quantity and found x ^ 10" '*-10'^ for typi- 
cal stellar populations (the only difference is that they calcu- 
lated the number of continuum photons just shortward of Ly a 
rather than shortward of 11.26 eV). Thus using our earlier 
definitions, we have 



ch 
47r 



Xnh{z) f*N^b A/coii, 



(17) 



where nh{z) is the mean baryon density and A/con = /coii(z) — 
/coii(-Zmax)- Neglecting rccombinations, wc Can then Write x; = 
^He/*/escA^7/j/coih SO the ionizcd fraction x„, when we reach 
the photodissociation threshold y„, 21 is 



;0.03 ( — 



J n 



.21 



0.1 



A/coll 



l+z 
10 



(18) 



Even with /esc — 1, X = 1' ^"d /cou ~ 3A/coii, the ionized 
fraction is at most several percent when we reach the pho- 
todissociation threshold. This estimate agrees with the de- 
tailed models of Haiman et al. (2000) (their Figure 7). Thus 
although this mechanism can cause x, to turn over, it can only 
do so near the beginning of reionization and the amplitude of 
the change will be small. 

Finally, we note that in reality the chemistry of H2 is suf- 
ficiently complicated that the detailed consequences of radia- 
tive feedback remain unclear. Initially it was thought that X- 
rays could catal yze its forma tion by increasin g the free elec- 
tron fraction (Haiman etal. 2000); however. lOh & HaimanI 
(|2p03) showed that the heating that inevitably accompanies 
X-rays impedes collapse, ensuring a net suppression. Others 
have argued for more c omplicated positive feedback rn echa- 
nisms near H II regions (iRicotti et al.ll200l ICenll2003ah . We 
will ignore all of these possibilities simply because positive 
feedback obviously cannot cause double reionization. More- 
over, the "self-regulation" process will not necessarily com- 
pletely halt the formation of low-mass halos. If small halos 
are responsible for the background, they will reach an equi- 
librium in which the UV background prevents the smallest 
halos from forming stars but allows more massive and better 
shielded halos to cool. Again, such a scenario tends to wash 
out double reionization. 

5. THE TOPOLOGY OF H II REGIONS 

One promising observable of the reionization process 
is the size distribution of H II regions and its evolu- 
tion with time (FZH04; Wyithe & Loeb 2004c), which 
can be probed through 21 cm tomo graphy (Furlanetto et al. 
|2p04c) or through Lya absorption (Furlanetto et al. 2004a; 
IWvithe & Loeb 2004b). By modulating the ionizing effi- 
ciency through local feedback mechanisms, the schemes we 
have described can in principle have substantial effects on 

^ We will assume for concreteness that all of the ionizing photons are con- 
tained in the frequency decade closest to the Lyman limit, which is a good 
approximation for stellar spectra. Even if the ionizing sources have hard 
spectra similar to quasars, our conclusion would not change. 



the size distribution of ionized bubbles. Here we will briefly 
discuss those consequences. Furlanetto et al.l (l2()04c) made 
some progress in this direction by considering simple scenar- 
ios with multiple generations of sources. We will improve on 
that treatment by examining the transition between star for- 
mation modes in more detail. 

We take the model of FZH04 as a starting point. We assume 
as above that the number of ionizing photons produced in a re- 
gion is proportional to the collapse fraction within the region. 
We denote the proportionality constant by C,. Note that this 
differs from the definition in ^2.4l in two ways: it is cumula- 
tive over the integrated star formation history and it includes 
a correction for past recombinations. A region can ionize it- 
self if C/coii > 1 ■ FZH04 showed how to transform this sim- 
ple condition into a size distribution by including the implicit 
scale dependence of /coii and ionizations from neighboring re- 
gions. We rewrite the ionization constraint as a condition on 
the density, 

5,n > <5.(m,z) = S,{z) - V2/:(C)[aL - a\m)]'^\ (19) 

where K{Q = erf^'(l — C~')- Equation il9i is then used 
as an absorbing barrier within t he ex cursion set formalism 
(Bond et al. 1991; Lacey & Co Ol993h to construct the mass 
function of ionized bubbles. The crucial point is that the 
shape of the barrier determines the distribution of bubble 
sizes. The a bove barrier is nearly linear in cr^, which implies 
jShethll99a) a well-defined characteristic size for the H II re- 
gions. Thus, the shape of the barrier is fixed by the condition 
/coll =constant, while the normalization is fixed by the ioniz- 
ing efficiency. 

We may now incorporate feedback by considering two sets 
of sources with different ionizing efficiencies and/or differ- 
ent galaxy mass thresholds. In the latter case (i.e., photoion- 
ization heating), the qualitative effects are easy to guess. A 
high-density region forms halos with Tvir ~ Tc which ionize 
the neighborhood. Without feedback, the associated bubble 
would grow through two processes: (/) it would merge with 
other H II regions, and (//) halos inside it would continue to 
accrete rapidly, producing more ionizing photons. However, 
if 7), ^ Tc, structure formation is suppressed in this region, the 
"internal" ionizations diminish, and the bubble can grow only 
by merging with its neighbors. Thus we would expect H II re- 
gions to be smaller and more numerous, because those around 
high-density regions grow more slowly and force voids to be 
ionized by the relatively rare galaxies embedded inside. At the 
same t ime, the high-density d ormant regions begin to recom- 
bine. Furlanetto et al.l i2004d) argued that this process would 
not stop until the entire universe had been ionized, though not 
necessarily simultaneously (see so eventually the pattern 
induced by recombinations (and exaggerated by the time de- 
lay between reionization in different regions) would dominate 
the signal. Unfortunately, because the photoheating depends 
directly on the ionized fraction, it cannot be incorporated into 
the FZH04 model. 

The FZH04 model can accommodate a transition due to 
metal enrichment more easily, because the ionizing efficiency 
does not depend on the local ionized fraction in this case. We 
will assume that there exists a function pii{6) characterizing 
the probability that any source in a region of overdensity S 
contains Pop II stars. Physically plausible models will have 
p[i > 0, where a prime denotes the derivative with respect 
to density, because structure formation is more advanced in 
dense regions. Given the ionizing efficiency of each star for- 
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mation mode, the condition for self-ionization is 

UniS) [(Cii - Cm)pii{S) + Cm] > 1 • (20) 
This modified barrier has two effects on the bubble size distri- 
bution. Unlike for a single type of ionizing sources, the left- 
hand side need not be a monotonically increasing function of 
the overdensity S. Its derivative vanishes when 



Pii 



-Pu 



Cii 



coll 



Cm — Cii 



(21) 



Physically, if ^ni S> Cii ^nd if high-density regions have en- 
riched themselves, then it is possible for these regions to re- 
main neutral even though regions with slightly smaller densi- 
ties (and thus incomplete enrichment) can self-ionize. Obvi- 
ously the increased recombination rate in dense regions will 
amplify this trend. 

As an example, we consider the choice pu = 1 — 
exp(— ?7/coii), which essentially reflects the wind model from 
^(ignoring /new)- Here the enrichment probability is a func- 
tion of density through its dependence on the collapse frac- 
tion; the exponential accounts for overlap of bubbles from 
randomly distributed sources within each region. Equation 
i2l\ then implies that the ionizing efficiency turns over at 
Vfcoii ~ 0.8 for (ill 3> Cii' ™d some regions sufficiently far 
above this threshold could remain neutral.^ Note that not ev- 
ery region with 77/coii > 0.8 would be neutral; above this point, 
the effective ( decreases with increasing S but not necessarily 
by enough to become neutral. Moreover, a dense region could 
still be ionized by its neighbors (Barkana & Loeb 2004). The 
turnover with density can mimic "outside-in" reionization in 
which voids are ionized before intergalactic sheets and fila- 
ments (where the sources are located). However, for this ef- 
fect to be significant we must have rifco\i{5 — 0) ^ 0.8; this 
only occurs when enrichment nears completion. Thus in this 
simple, wind-driven model we do not expect neutral dense re- 
gions to be common during reionization. 

Unfortunately, once the ionizing efficiency turns over, the 
excursion set formalism of FZH04 breaks down. This is be- 
cause the number of ionizing photons is no longer additive 
in a simple way, but rather depends on the degree of enrich- 
ment of substructure inside the region of interest. However, 
if r]fco\\{S = 0) remains reasonably small during reionization 
then the model can still be used. In this limit, we consider 
how the shape of the barrier - and hence the shape of the bub- 
ble size distribution - changes due to the metal enrichment 
pattern. In standard reionization scenarios, the shape is set 
by /coll ^constant (the solution to C/coii = !)■ With enrich- 
ment, equation (I20> shows that the shape is instead set by a 
condition on a combination of /coii and pn. Unfortunately, 
the dependence is weak for reasonable choices of pu. For ex- 
ample, with the simple wind model examined above, pu is 
purely a function of /cou and the shape is still determined by 
/coll =constant; it can be exactly reproduced by choosing an 
effective ^etf to match the total ionized fraction. 

Thus, to change the shape, pu must have a more compli- 
cated density dependence. Let us consider the simple param- 
eterization 



pii^l- exp 



1+A. 



l+z 



(22) 



where A is some constant that measures the steepness of the 
density dependence. This form has the limits pu 1 for 

* Note that f^^u actually depends on the scale of interest (as could pu)- We 
neglect this here, considering only those scales for which (T,^i„ 2> cr(m). 




Fig. 5. — Sample barriers from equation 1221 (solid curves) and for a single 
type of source (dashed curves) with the same limiting value at = 0. In the 
solid curves, C,i\\ = 44, Qi = 9, and A = 1. We show results for rj = 5, 15, 
and 25, from bottom to top. 



S ^ 00 and pii for 6 —00, and it concentrates enrich- 
ment in dense regions more than the simple wind model does. 
One reason this may happen is that /new//coii is smaller than 
the average in dense regions because they are further along 
in the structure formation process. Qualitatively, we expect 
the barrier to be somewhat steeper with this prescription, be- 
cause dense regions have extra enrichment and fewer ionizing 
photons. The solid curves in Figure |5] show some example 
barriers. In each case, Cm = 44, (u = 9, and A = 1. We 
show T] — 5, 15, and 25, from bottom to top. The dashed 
curves show the barrier for a single source type with the 
value at cr^ = fixed to match the solid curves. These have 
Ceff — 38, 24, and 1 1 for the three values of 77. Indeed the bar- 
riers are steeper, which implies fewer small bubbles because 
small, dense clusters are highly enriched. However, we find 
that the steepening is relatively weak, especially when one of 
the two star formation modes dominates (i.e, Ceff ~ Cm or Cn)- 
Of course, if enrichment is sufficiently extensive that equa- 
tion i2l\ is violated, the barrier approach breaks down and 
the bubble pattern could change significantly. 

We have considered other possible parameterizations of pu 
based around the simple wind expression; in the cases we ex- 
amined, the differences were not much larger than shown in 
Figure Is] The reason is that pu is simply used to interpolate 
between two different populations, each of which would fol- 
low the usual barrier In addition, the interpolation between 
them also depends on the collapse fraction. Thus it appears 
difficult to change the shape dramatically, although exotic sce- 
narios in which the collapse fraction played no role could have 
stronger effects. 

6. DISCUSSION 

We have examined a variety of feedback mechanisms that 
could potentially cause "double reionization," in which the 
globally-averaged ionized fraction x, decreases over a limited 
period of cosmic time. These mechanisms include metal en- 
richment, photoionization heating, and photodissociation of 
H2. We constructed simple models that nevertheless retained 
the crucial physics; as a result we were able to vary a rela- 
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lively small set of parameters to pin down the requirements 
for double reionization. In all cases, we found that double 
reionization requires extreme (though not impossible) param- 
eter choices. In particular, it requires a rapid drop in the ion- 
izing emissivity over a single recombination time (see equa- 
tion 1121 '). Because the recombination time varies spatially 
within the IGM, we found that it is not difficult to imagine 
relatively dense pockets recombining after feedback becomes 
strong (which could affect the topology of ionized gas; see 
but because these occupy only a small fraction of the 
universe, double reionization is difficult to arrange on large 
scales. 

For metal enrichment, we require three conditions: 
(0 Ciii/Cii ^ 1' (") ™ average clumping factor sig- 
nificantly larger than unity, and (///) winds that ex- 
pand much more rapidly than predicted by existing semi- 
analytic models ( Madau et al. 2001; Scannapieco et al J2002t 
iFurlanetto & Loeb 2003). We have argued that con- 
dition (//) is implausible for the majority of the IGM 
(jMiralda-Escude et al. 2000') and that (///) is similarly implau- 
sible, so we regard this mechanism as unlikely to cause dou- 
ble reionization. Existing work predicted double reionization 
only because of an artificial, instantaneous transition between 
star- forming modes ( Wvithe & Loeb 2003a; Cen 2003a' see 
also 'Haiman & Holder 2003). Our models are more similar 
to Scannapieco et al. (2003), who found a gradual disappear- 
ance of Pop 111 stars. In a hierarchical picture, the decline 
has two causes: accretion shifts to higher mass, pre-enriched 
halos even when the global metallicity is small and the slow 
expansion of galactic winds. The first depends on unknown 
halo merger rates, but our estimate appears to be conserva- 
tive. 

Double reionization from photoheating is similarly diffi- 
cult. In this case we require (;) TijTc ^ 1, (//) small halos 
(those with Tvii < 7),) contribute a large fraction of the ion- 
izing photons, and (///) new ionizing sources form in regions 
that have not yet been ionized, so that their photons are not 
"wasted" on recombinations (because in that case they do not 
heat the IGM). The second condition implies that the star for- 
mation efficiency does not drop too rapidly with decreasing 
halo mass. In particular, if cx m^/^ persists at high red- 
shifts (as has been observed in nearby galaxies; Dekel & Woo 
l2()03tlKauffmann et alJl200 3l). then small halos do not domi- 
nate the ionizing photon budget and double reionization can- 
not occur Quasars also help to wash out double reionization 
if the black hole mass scales steeply with host halo mass, as 
in Wvithe & Loeb (2003c). We argued that the photoheating 
occurs more rapidly than ionization, because photons are pref- 
erentially emitted in regions that have remained neutral to that 
point, while recombinations are confined to ionized regions. 
We made an extreme assumption in order to maximize the 
likelihood of double reionization; this picture can be tested 
with self-consistent photoheating feedback in simulations. If 
Xh ~ X; is more accurate, double reionization is impossible to 
achieve through photoheating. Our fiducial choice for the de- 
gree of suppression, Th/Tc = 25 is on the threshold of double 
reionization; if photoheating is less efficient at supp ressing 
galaxy formation near the time of reionization (.Diikstra et al.l 
|2Q04a), then this mechanism also cannot cause double reion- 
ization. Note also that, if the first sources produce X-ray pho- 
tons (e.g., Ricotti et al. 2004; but see Dijkstra et al. 2004b), 
the IGM will be heated gradually and (nea rly) uniformly, 
again suppressing accretion onto small halos dOh & HaimanI 
120031) . This would tend to smooth out the transition from Tc to 



Ti„ decreasing the likelihood of double reionization. Clearly, 
double reionization through either metal enrichment or photo- 
heating requires the confluence of a number of unlikely pos- 
sibilities. 

Finally, H2 photodissociation has a similar phenomenolog- 
ical effect to photoheating, increasing the minimum mass for 
star formation (albeit from a much lower base level). How- 
ever, we showed that the feedback threshold must be reached 
when the ionized fraction is still small, so the recombina- 
tion phase (even if it occurred) would be difficult to observe. 
iHaiman et a l. (2000) also reached this conclusion through 
more detailed modeling. 

We have not examined any combinations of these different 
mechanisms. Together they could increase the contrast be- 
tween the initial and final ionizing efficiency, which would 
make double reionization somewhat easier to achieve. For 
example, photoionization heating could extend reionization 
over a long interval. This would increase the time over which 
winds can expand and hence move the Pop III/II transition 
closer to reionization, increasing the effective amount of sup- 
pression in ionized regions. Another possibility is that star 
formation through H2 cooling at extremely high redshifts 
could begin to spread metals throughout the IGM. Once pho- 
todissociation dominates, there would be a long pause until 
reionization could continue. Again, the winds would have 
extra time to expand. However, we see no reason for any 
such scenario to increase qualitatively the likelihood of dou- 
ble reionization. 

We have focused here on particular ionization histories in 
which X, turns over We found that such scenarios are diffi- 
cult to achieve. But we stress that extended reionization is not 
nearly so difficult; any form of feedback will help to prolong 
the reionization era and to relieve the apparent tensi on be- 
tween t he z ^ 6 quasar da ta (Fan et al. 2002; Wvithe & Loebl 
"2004a* 'Mesinger & Haiman 2004) and the optical depth to 
electron scattering measured by WMAP (Kogut et al. 2003). 
Indeed, a number of our models in Figure s|3] and |4] show that 
feedback can significantly extend the era over which x,- ~ 0.5. 
Moreover, we have found that a turnover in the emissivity is 
also not difficult to achieve; many of our feedback scenarios 
yield this kind of behavior, at least if low-mass halos have 
high star formation efficiencies. Unfortunately, the emissivity 
evolution is harder to measure. 

Finally, we must also stress that we have studied double 
reionization in a global sense. Although we have found that 
x, only turns over in exceptional circumstances, the reioniza- 
tion history of any particular region or line of sight need not 
be monotonic. For example, if photoheating suppression is 
important, a volume ionized early on will experience a sharp 
decrease in its emissivity and begin to recombine. Only later 
will halos grow sufficiently massive to ionize it again. Our 
models show that, while this region is recombining, others 
are also being ionized, and the latter generally dominate the 
global evolution. The most promising way to probe reion- 
ization in different physical volumes is 21 cm tomography 
(Madau etal. 1997; Zaldarriaga et al. 2004; Furlanetto et al] 
|2()04c). We have argued that these feedback mechanisms 
should induce signatures in the bubble size distribution, al- 
though metal enrichment appears to have little effect except 
in dense, highly-enriched regions. Because of the complexity 
of this problem, numerical simulations would provide the best 
method to examine changes in the bubble pattern in greater 
detail. 
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APPENDIX 

AN EXTENDED PRESS-SCHECHTER ESTIMATE FOR THE HALO FORMATION RATE 

In ^2.11 we estimated d/new/dz in a simple fashion. Here we point out that self-consistent analytic calculations for this quantity 
do not yet exist. Given the comoving number density of halos per unit mass as a function of redshift, n{m,z), and the rate at which 
these halos merge per unit volume, «(mi,z) «(m2,z) Q{mi,m2,z), we can compute the r ate at which halos merge to form new 
objects with masses abov e OTrnin- Unfortunately, as recently em phasized by Benson et al.l J2004h . the merger rates typically used 
in cosmological studies ( Bond et alJl99lHLacev & Colel 19931) are not well-defined. The problem is most easily seen by noting 
that Q{ini,m2) ^ Q{m2,nii) in the extended Press-Schechter (EPS) formalism. The fundamental difficulty is that this formalism 
does not assign points in space to unique halos: the smoothing procedure implicit in these approaches requires that neighboring 
points be identified with halos of different masses, even though they must physically be part of the same object. Merger rates, 
which by definition require uniquely specified halos, are thus ill-defined quantities. As a consequence, merger tree algorithms 
based on the EPS formalism either do not reproduce the proper mass function or do not conserve mass during the individual time 
steps, although ad hoc procedures to alleviate these problems have been developed (e.g., Somerville & Kolatt 1999). 

Nevertheless, it is useful to compute d/new/dz with the EPS formalism in order to investigate whether the slow evolution in 
equation (|3 is an artifact of our approximation. A co mparison to the EPS estimate is motivated by its succ es in areas such as 
semi-analytic galaxy formation (e.g., Kauffmann & Whitdll993tlCole et alJll99'4llSomerville & Kolattll 1999ft and modeling the 
quasar luminosity function (e.g., Wvithe & Loeb 2002). 
For a halo of mass m at redshift zi, the EPS formalism provides the fraction of matter that was in objects with masses smaller 
at some earlier time Zh iLacev & Colel 19931) and hence the fraction of a halo's mass that was accreted recently: 



^'(< m^m,Zh\m,zi) =erf 



Scizh) - 6cizi) 



(Al) 



The relevant timescale is the dynamical time within a galaxy (which gives approximately the time over which star formation 
occurs), fdyn ^ {Gpy^/^ ^ (A,.pcG)"'/2 ^/4/277^j /~'(^), where p,. = 3//2/87rG is the critical density, H{z) is the Hubble 
parameter at redshift z, and A,, is the virial overdensity (iBarkana & Loeb 2001 ) . If we then assume that accretion onto halos with 
m < 2mniin can be included in the newly-formed halo component, we can estimate the fraction of collapsing gas that accretes 
directly onto existing halos via 

/2m,„i„ dm m n{m) F {< m^in,Zh\m,zi) 
/old — noo , r^7r~ j Tj (A2) 

j^^^^ dmmn[m)F(< m^^,Zh\m,Zi) 

so that d/new/ dz w ( 1 — /oid)d/coii / dz. The denominator is necessary to normalize the total mass accretion rate to its actual value; 
a direct calculation with the EPS merger rates of the total collapse rate does not reproduce the true d/cou/dz. One advantage of 
this formulation is that it includes the duration over which a new halo must remain isolated (unlike in equation Q). In other 
words, a halo that merges into a massive object immediately after passing above m^in will most likely form Pop II, rather than 
Pop III, stars. We have used the dynamical time to fix Zh- Varying the time offset has only a small effect on the results: increasing 
or decreasing f^yn by a factor of two changes d/new/dz by < 10% and leaves the shape unaffected. Another approach is to use 
the explicit expression for Q{mi,m2,z) in the EPS formalism to compute d/new/ dz di rectly (but without the time offset). This 
procedure (with symmetrized EPS merger rates) yields results similar to equation \A2\ . 

The resulting formation rate of new halos is shown in the bottom panel of Figure^for the same three mass thresholds as in the 
top panel. Interestingly, in this case d/new/dz evolves even more slowly than our standard calculation predicts: it falls below unity 
earlier and remains substantial much later. We thus conclude that d/new/dz, while unknown in detail, has the important property 
that it evolves slowly and smoothly with redshift, vanishing only at z < 6 for the mass thresholds of interest. 
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